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Abstract 

We study the phase diagram of the 4d compact U(l) gauge theory as a function 
of the number of Euclidean time slices. We use the helicity modulus 1 as order 
parameter to probe the phase transitions. The order of the transition along the 
phase boundaries is studied and the possibility of a continuum limit is discussed. 
We present new, strong evidence that the T = bulk phase transition is first-order. 
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1 Introduction 



The way in which continuum QFT incorporates the definition of 'temperature' has 
many relations with the notion of dimensional reduction: one crucial ingredient is 
in fact the compactification of the temporal dimension, whose length is tuned as a 
free parameter according to 



where T is the temperature and | L t | is the (dimensionful) time extension. On the 
lattice Eq.(JTJ) reads 



where a is the lattice spacing and now L t is a pure integer number; in this case the 
notion of finite temperature requires the possibility of taking the limit a — > 0, and 
then to tune L t accordingly to keep T fixed. 

In this paper we consider the 4d compact U(l) gauge theory on a lattice at fixed 
L t . We study its phase diagram and pay particular attention to the points at which 
a continuum limit might be extracted, namely the phase boundaries. This is an 
interesting generalization of the analogous problem at 'zero temperature' (L t = oo), 
still open from the theoretical point of view, where only recently convincing numer- 
ical evidence of the first order nature of the transition was obtained QjJ E] . Further 
investigations seem therefore well motivated. 

The paper is organized as follows. In Sec. El we describe the model and fix 
the formalism. In Sec. H3 we introduce the definition of the order parameter and 
review some of its properties. In Sec. 0]we study the phase diagram of the model. 
Conclusions follow. 



We consider a four dimensional lattice ■ L t with periodic boundary conditions. 
We associate to each link U^r) (r = (x, y, z, t)) connecting two neighboring sites an 
element of the group U(l), namely e %e ^ r \ where 0^{r) = aeA^{r) is proportional to 
the vector potential A M (r) and e is the bare gauge coupling. The system is governed 
by the Wilson action 
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Figure 1: Cartoon of a temporal (s, t) plaquette in the case L t = 1. The two spatial 
links cancel each other, leading to a decoupling of the temporal links. 

where {9p{x)) llv = 9^(x) + 9 v (x + fi) — Q^ix + v) — 9 v (x) is the so-called plaquette 
angle, and j3 — -y. 

The phase structure of this model is well known at zero temperature (L t = 
L s = L, L — > oo): it has a strong coupling confining phase for (3 < j3 c , and a 
weak coupling Coulomb phase for j3 > f3 c . There is strong numerical evidence that 
the phase transition (at (3 C = 1.0111331(15)) is of first order [T]|2] . Thus, in the 
extended phase diagram (fl, 1/L t ), we have a complete picture of the axis l/L t = 0. 

There is another limiting case which is easy to understand, namely L t = 1: 
consider one temporal plaquette, represented in the cartoon of Fig. [T] In the case 
Lt — 1 the temporal link is already a gauge invariant quantity, a Polyakov loop 
itself. Moreover, from Fig. ^ it is clear that around a temporal plaquette the two 
spatial links are actually the same link considered twice in opposite orientations, 
therefore they exactly cancel each other. The temporal plaquette angle is therefore 

e s (x) + e t ( x + §)- e s (x) - e t (x) = e t (x + §)- e t (x) (4) 

The contribution of the temporal plaquettes to the action is thus 

Stomp. = ~I3 ^2 cos 9 p = -P/, cos(9 t (x + s) - 9 t {x)) (5) 

-Ptemp. X 

and it is completely decoupled from the spatial plaquettes. Furthermore, Eq.© is 
the action of the 3d XY model. One then realizes that in the case L t = 1 the 
partition function of the model factorizes into 

ZLt=l — %3d XY ' Z 3d LGT (6) 

namely we decompose the full partition function into that of the 3d XY model and 
that of the 3d U(l) lattice gauge theory. We know that the former undergoes a 
phase transition of second order at /? = 0.45420(2) [3], while the latter is always in 
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Figure 2: Sketch of the phase diagram of the model, collecting the information 
about the two limiting situations L t = 1 and L t = oo. 

the confined phase. We summarize the information gained so far about the phase 
diagram in Fig. |21 

In order to explore the phase structure of the model for a generic value of L t we 
need an order parameter; we introduce it in the next section. 

3 The order parameter 

To distinguish between confining and Coulomb behavior, we could monitor the ex- 
pectation values of large Wilson loops, or the correlation of two Polyakov loops. 
However, we can also construct an observable which directly probes the large- 
distance properties of the system, and which has proven itself to be an efficient 
tool for analyzing a weak first-order transition Jlj: the helicity modulus. 

To introduce the helicity modulus for this theory, we first review the notion of re- 
sponse function to an external static field. Consider our four-dimensional lattice and 
choose an orientation (//, u); if we call L p , L a the two other (orthogonal) directions, 
a stack of L p ■ L a parallel (/x, v) planes is thus defined. Let us modify the partition 
function of the system by imposing an external flux $ through this orientation, as 
shown in Fig. El the flux is homogeneously spread over each of the parallel planes, 
meaning that $ P = &/L p L u is assigned to each plaquette with orientation (/i, u), 
as written in the following equation 

S($) = -P cos(fl P + z -^)-/3 J2 cosdp W 

P&{^i,y orient.) " U P^(p J ,u orient.) 
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Figure 3: Graphical representation of the static external flux imposed through the 
(}Ji,v) orientation. 

One may wonder if this specific choice of spreading the external flux plays a role 
in our construction, but the answer is no, since in this Abelian theory a suitable 
change of variables flj can vary at will the distribution of the flux, leaving the 
partition function of the system unchanged. Moreover, one can consider another 
point of view, and show that the external flux can actually be imposed also via a 
change of the boundary conditions from pure periodic to twisted: this is very easily 
accomplished if we impose, e.g., that each link U^ir) in a particular // direction is 
identified with its image at distance L u up to a phase factor e % ^l Lii : 

9^x + L v ) = d^x) + ^- (8) 

A* 

Let us now consider the free energy of the system 

= -logZ($) = - log J dUe- s{u ^ (9) 

and study its behavior in the two phases. In the confining phase the correlation 
length £ is finite (a natural scale in this phase is provided by the string tension 
a = l/£ 2 ), and the effect of the external flux, i.e., of the boundary conditions, 
decreases as e~ Lv ^. Therefore, in the thermodynamic limit, we obtain the result 

= const. V $ (10) 

The situation is different in the Coulomb phase, where the correlation length of 
the system is infinite (the Coulomb potential does not fix any scale in the system) 
and therefore the free energy of the system has a non-trivial dependence on <3>. In 
order to understand quantitatively this dependence, consider Eq.(JZJ) in the classical 
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limit (3 — > oo, where all fluctuations are suppressed and the action of Eq.(J7j) reduces 
to the following 

^classic ($) = ~P }^ COS ( -^-^ )+ COnst. (11) 

P£(/i,v orient.) " 

namely, only the contribution of the twisted plaquettes is present. Consider the 
thermodynamic limit at finite temperature, i.e., keeping one size fixed and letting 
the others diverge. The product ■ L u diverges, whatever orientation we choose, 
therefore it is possible to expand Eq.([ll|) around its argument, and obtain 

Fclassic($) " Fclassic($ = 0) = (12) 

since the sum is over V = L^L v L p L a identical contributions. Some comments about 
this last equation are necessary: first of all, we showed in P that for a finite (3, where 
fluctuations are present, Eq.()12p is only slightly modified, namely 

infinite " ^nitc = 0) = M^Ml ^ 

i.e. the bare coupling is replaced by a renormalized coupling, and this is all we need 
to effectively take into account the fluctuations of the system. Another remark is 
that Eq.(J7|) is 2ir periodic in $, while Eq. (jl3j) is not. The corresponding modification 
of Eq. (|13|) is exponentially small in the flux $ for small $ pQ , and plays no role in 
what follows. 



Reconsider now Eq. (jl3|) . and define 

HP) = |*=o (14) 

This quantity is called helicity modulus in the literature, and it is a measure of 
the curvature of the free energy profile around $ = 0. From our previous discussion 
about the behavior of F(<&) in the two phases, it follows that h(P) is an order 
parameter for our theory. If one performs explicitly the double derivative w.r.t. $, 
one finds 

&Q3) = - r L-{< P^s6 P }-(( P^ep) 2 )} (15) 

^ U P£(fi,u orient.) P(L{[l,i> orient.) 

directly amenable to numerical simulations. 

It is important to observe that in our finite temperature system we expect a 
different response of the system to fluxes imposed through spatial or temporal ori- 
entation. Consequently the helicity modulus is expected to show different behavior. 
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Figure 4: Conjectured phase diagram of the model. The phase boundary on the left 
separates a confining phase from a partially Coulomb phase, in which static electric 
charges propagating in the temporal direction interact Coulombically while magnetic 
charges propagating in a spatial direction are confined. The phase boundary on the 
right separates this partially Coulomb phase from a pure Coulomb phase. The 
dotted line indicates the possibility that the transition remain second order for a 
range 1 < L t < L t of temporal sizes. 

In the following we will carefully distinguish the two situations, calling spatial (resp. 
temporal) helicity modulus the response to a flux through spatial (resp. temporal) 
planes. 

4 Phase structure of the model 

Let us go back to Fig. |2] and try to make an educated guess about the phase structure 
of the model. First of all, we expect a phase boundary connecting the first order 
transition point at L t = oo with the second order point at L t = 1 (see Fig. HJ). 
On the left of this curve we are in a confining phase; on the right we expect a 
different behavior of temporal and spatial helicity moduli. A natural possibility is 
that the T = transition between the confined and Coulomb phases splits for the two 
orientations at finite temperature. As a consequence we would have an intermediate 
phase in which temporal Wilson loops obey perimeter law, while spatial ones would 
still show confinement (area law), just like in the high temperature phase of Yang- 
Mills theories. 

Another important issue is the order of the phase transition along the phase 
boundary on the left: we know that in the two limiting cases (L t = 1 and L t = oo) 
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the order of the transition is different. Two scenarios are possible: either the transi- 
tion remains more and more weakly first order, and turns to second order only when 
Lt = 1, or there is some special value L t > 1 for which the transition turns already 
to second order. This second scenario is represented in Fig. EJby the dashed line. 

Let us now consider and analyze separately the two conjectured phase bound- 
aries, to understand the structure of the model as the thermodynamic limit is ap- 
proached. 

4.1 Transition to the Coulomb phase 

In this section we consider the phase boundary on the right in Fig. the one 
connecting the temporal Coulomb phase with the pure Coulomb phase. To explore 
this transition we want to observe a change of behavior (from area to perimeter 
law) in the spatial loops, therefore we have to probe the response of the system 
to external fluxes through spatial planes. Equivalently, we will consider the spatial 
helicity modulus. 

In Fig.EJwe show the numerical results for the spatial helicity modulus, measured 
on a set of lattices with temporal size L t = 2 and different spatial sizes L s . What is 
apparent is that the transition point does not seem to converge at all to some fixed 
value /3 C , but moves to higher and higher values with the spatial lattice size. If this 
is the case, this transition disappears from the phase diagram in the thermodynamic 
limit. 

A first rough way to understand the behavior of the system in this case is given 
by Eq.([T3|). Let us rewrite it here, focusing our attention on the geometrical factor 

^finite - ^finite = «) = M^^L^L ( 1 6) 

where we have highlighted the role of the spatial and temporal directions in this 
case. At any fixed value of $ and at any finite (3 we get therefore 

- F(0) ~ ^ -> as L s -> oo (17) 

meaning that the flux free energy always vanishes in the thermodynamic limit; this 
prediction is nicely confirmed in Fig- El where F($>) is measured from the flux distri- 
bution z/($) = e~ F ^ for different spatial extensions L s . As a consequence, no phase 
transition is present between the mentioned phases for any finite value of L t , but 
just the interplay of the coupling and of the spatial lattice size mimics a transition. 

An alternative way to understand this finite size effect is the following. At any 
given value of f3 there is a finite density of magnetic monopole currents. In particular, 
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Figure 5: Spatial helicity modulus as a function of (3, measured for a temporal 
size L t = 2 and for different L s . The value of (3 at the apparent transition point 
grows with the spatial size L s . In the thermodynamic limit this phase boundary 
disappears. 

a wrapping (non-contractible) time-like current loop, far apart from its companion 
anti-loop, can disorder all the spatial planes in the lattice, or, equivalently, a spatial 
Wilson loop of arbitrary size; this, in turn, produces confinement [Tj. Call p(f3) the 
density of currents at a given value of /3; we know j^j that 



where c is a constant, for (3 1. Let n(l)dl be the distribution function indicating 
the fraction of monopole loops of length between / and / + dl. The transition occurs 
approximately when each configuration contains a pair of non-contractible monopole 
loops, namely when 



where the integration domain I > 2L t represents at least a necessary condition to 
have wrapping loops; on average we expect it to be also sufficient. The last piece of 
information we need is 



p{(3) ~ e 



-c/3 





(19) 



dn(l) 
dl 



-si 



(20) 



8 



8 3. 2 . 

12^*2 x 

16^*2 + 
24,*2 

32 J *2 ° 



+ ^=hu +^ + 

" + X + X x X + x 4 " 

■ ■ • 

□ □ 
□ 1? # □ 



-4jc 



-2ti 



2tu 



4rt 



Figure 6: Flux free energy for lattices with temporal size L t = 2 and increasing 
spatial size, at (3 = 1.4. As L s — >• oo the flux free energy tends to zero, as expected 
from Eq. (jl6|) . 

motivated by the fact that the free energy cost of a current is proportional to its 
length / (notice that this is true only in the Coulomb phase, where the coefficient 
5 > 0). Let us put together all the ingredients and rewrite Eq. (fn?j) 

poo p-2&Lt 

L\e~^ \ e~ 5l dl = L\e-^ e -— > 1 (21) 

which contains all the information we need. At fixed volume, the inequality is satis- 
fied only for sufficiently low (3\ when (3 is large enough we are always in an ordered 
phase, as indicated by our numerical finding (spatial helicity modulus different from 
zero). In the limit L s — > oo at fixed L t , this equation is always satisfied, indicating 
that the system is always disordered. It is only in the limit L t — > oo (because of 
the term e~ 2SLt ) that we can obtain an ordered phase. Interestingly, this equation 
predicts a pseudocritical coupling 

f3c~ log L s (22) 
contrary to the simpler argument in Eq. (|17jl which gives 

Pc ~ L s (23) 

The first argument is completely classical, therefore we expect that Eq. (j22j) describes 
the physics better. 
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Figure 7: Temporal helicity modulus as a function of /3, measured on systems of 
temporal size L t = 4 and different spatial sizes L s . 

4.2 Transition to the confining phase 

Let us now turn our attention to the left phase boundary in Fig. HI and to the 
temporal helicity modulus. Eq.(fT5j) gives 

ignite M $ ) " ^[finite = 0) = ^l^LsLs ^ LgW) (24) 

therefore in the thermodynamic limit, at the transition point, the signal is infinitely 
enhanced; the transition point is determined (up to finite size effects) by the be- 
havior of /3r(/3), which is related to the free energy of a monopole current, an order 
parameter by itself. In Fig. 0we show one set of results for L t = 4 and three spatial 
volumes (the data are reweighted a la Ferrenberg-Swendsen |.6j). From a finite size 
scaling analysis of such data we can extract the numerical value of (3 C as a function 
of the time extension L t , and the order of the transition. For small values of L t we 
want to decide if the transition is first or second order. Our strategy is to assume 
the second order scenario (the universality class is then that of the 3d XY model, as 
for L t = 1) and then check the quality of the data collapse. We thus fix v = 0.6723 
(the value of the 3d XY model [3]) and rescale the (3 axis according to 

/3-+L 1 J v (p-p e )~fe) 1 ''' (25) 

where we have to tune f3 c only. The results of such analysis for the lattices with 
L t — 2, 4 are presented in Fig. |H| Comparing the excellent collapse of the curves 
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the temporal helicity modulus for lattices 
= 4 (right). The scaling window becomes 
the approach to a first order transition. 



in the case L t = 2 with the mediocre collapse in the case L t = 4, one gets a first 
qualitative indication of a second order transition turning into a first order one. 
Moreover, one can speculate that for small values of L t > 1 the transition still is 
second order, as indicated in the cartoon of the phase diagram of Fig. 0] Since 
this issue is relevant to the possibility of extracting a continuum limit, we further 
investigated this problem directly by measuring the latent heat of the transition 
as a function of L t . We took the direct approach of monitoring the Monte-Carlo 
history of the plaquette, starting from hot and cold (disordered and ordered) initial 
configurations. When a metastability is observed, we measure the distance between 
the two plateau values. We repeated the measurements for different volumes in 
order to keep the finite volumes effects under control. The plaquette histories are 
displayed in Fig. and the estimated latent heats AE(L t ) in Fig. [T0J Table I 
summarizes our measurements of the critical parameters. As far as we know there 
are no theoretical predictions about the functional form that we observe for AE(L t ). 
The data anyway seem to suggest a scenario in which the latent heat drops to zero 
already for L t ~ 4, at (3 c {L t ) ~ 1. 



We cannot exclude a small but non-zero latent heat of course. However, Eq. (J24j) 
suggests that the transition turns from first- to second-order when (3r(/3 c ) vanishes 
as the transition point f3 c is approached from above, since then the magnitude of 
the jump in the helicity modulus at the phase transition also vanishes. The rapid 
decrease of (3r with (3 observed at zero temperature in [T] makes it unlikely that 
/3r would remain non-zero all the way down to f3 c (L t — 1) « 0.45420. Thus, we 
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Figure 9: Time history of the plaquette, starting from an ordered (upper curves) and 
disordered (lower curves) initial state. The latent heat is defined as the difference of 
the values of the two plateaus. From left to right the cases L t = 8, 6, 4 are studied. 
For L t = 8, 6, where the latent heat is certainly different from zero, we checked that 
the thermodynamic limit had been reached by comparing different spatial volumes. 
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Table 1: Latent heat AE and transition coupling j3 c versus L t . AE decreases very 
quickly for L t < 6, opening the possibility of a second order phase transition for 
Lt <4. 



consider the scenario where the transition is second-order for a range of temporal 
sizes [1, L t ] as a likely alternative. 

This opens the intriguing possibility of a continuum limit, obtained by approach- 
ing any point along such a second-order line. Note however that the number L t of 
time-slices cannot exceed L t , so that one cannot strictly speaking talk about the 
£7(1) system being at "finite temperature". Rather, we have a system of L t < L t 
coupled 3d— systems. Moreover, as one approaches the critical point, only the tem- 
poral correlation length £ t /a diverges in lattice units, whereas the spatial correlation 
length £ s /a remains finite. Thus, keeping the "physical" correlation length £ t con- 
stant would imply not only a — > 0, but also £ s — > 0, so that the spatial planes would 
be completely disordered. The physical relevance of this situation seems marginal. 

The data in Table I can also be used for another kind of analysis, namely to study 
the dependence of (3 C on L t as L t — > oo. The coupling at the transition should in 
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Figure 10: Latent heat AE of the transition versus 1/L t . The displayed data 
correspond to L t = oo, then L t = 8, 6, 1. The cases L t = 4 and 2 are not displayed 
because the question whether the latent heat is already zero cannot be resolved with 
our numerical measurements. The dashed line corresponds to the simplest, linear 
ansatz. 

fact approach its asymptotic value (3 C {1/L t = 0) with typical finite-size corrections: 

P c (L t ) « P c (l/L t = 0) + (26) 
where c is a constant, and v = -? = ~ for a first-order transition. This is indeed the 

a 4 

case, as shown in Fig. El where the straight line corresponds to \jL\ corrections. 
This represents a qualitatively new and strong evidence of the first order nature of 
the transition in the zero temperature limit. 

5 Conclusions 

In this paper we have investigated the phase diagram of the 4d U(l) lattice gauge 
theory in the extended plane (J3, 1/Lt), where L t is the number of time-slices of the 
system. A summary of our findings is provided by Fig. IT2l where the phase diagram 
of the model is shown. Only one phase boundary is present, whose position we 
have determined. It divides a confining phase on the left and a temporal Coulomb 
phase (with broken global U(l) symmetry) on the right for any finite L t ; the pure 
Coulomb phase exists only for L t = oo and j3 > (3 C = 1.0111331(15). This phase 
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Figure 11: Shift of j3 c as a function of the Euclidean time extent L t . The scaling 
ansatz shown is —^7, with c = 2.18 and v — 1/4. An exponent v — 1/dis the typical 

signature of a first-order phase transition. 



diagram is remarkably similar to what we find in Yang-Mills theories, where spatial 
confinement persists at all temperatures. 

The order parameter we used, the helicity modulus, reduces naturally, via di- 
mensional reduction, to the ordinary helicity modulus known in condensed matter. 
We can therefore say that our order parameter is a generalization for gauge theories 
of the already known quantity. 

The transition for any finite value of L t is first order, except perhaps for L t < 4. 
In this case the transition becomes so weak that we cannot determine its order, but 
consider likely the possibility of a second order phase transition. This question could 
be studied by more powerful numerical methods. For instance, one could implement 
spatial C-periodic boundary conditions, which reduce the global U(l) symmetry to 
Z2, then measure the height of the free energy barrier between the two vacua as a 
function of the system size. 

However, even if the transition is confirmed to be second order for L t < 4, the 
relevance of this result to a continuum limit of the finite temperature U(l) theory 
is marginal: in order to introduce the notion of temperature 



T 



L t a 



(27) 



one must tune L t as a function of a — > in order to maintain a finite T. If, as in our 
case, Lt cannot be larger than a certain L t , the continuum limit cannot be taken, 
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Figure 12: Phase diagram of the model. The phase boundary is obtained from 
a finite-size scaling analysis of the temporal helicity modulus. Three phases are 
indicated: a confining phase on the left of the phase boundary, a temporal Coulomb 
phase on the right of the boundary at finite L t , and the pure Coulomb phase on the 
right of the boundary and for L t = oo. 

and we should not really talk about "finite temperature" . 

Finally, we would like to mention the old paper [7] , where the authors considered 
the Georgi-Glashow model with gauge group SU (2) at finite temperature; this model 
reduces to our U(l) model in the limit of complete gauge symmetry breaking (where 
U(l) is the symmetry remnant of the broken SU(2)). The conjecture of the authors 
about the phase diagram of our model, contained in their Fig. 1, is completely 
consistent with our findings. 

An even older, correct, phase diagram can be found in The persistence of 
the spatial area law at finite temperature was first proven in 
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